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Abstract 

We discuss the spectrum of screening masses in a pure gauge the- 
ory near the deconfinement temperature from the point of view of the 
dimensionally reduced model describing the spontaneous breaking of 
the center symmetry. Universality arguments can be used to predict 
the values of the mass ratios in the scaling region of the deconfined 
phase when the transition is of second order. One such prediction is 
that the scalar sector of the screening spectrum in SU(2) pure gauge 
theory contains a bound state of the fundamental excitation, corre- 
sponding through universality to the bound state found in the 3D 
Ising model and </> 4 theory in the broken symmetry phase. A Monte 
Carlo evaluation of the screening masses in the gauge theory confirms 
the validity of the prediction. We briefly discuss the possibility of us- 
ing similar arguments for first order deconfinement transitions, and in 
particular for the physically relevant case of SU(3). 
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1 Universality prediction for the screening 
masses 



QCD predicts the existence of a deconfinement temperature above which the 
fundamental degrees of freedom cease to be confined in color singlets and 
can propagate in isolation. This prediction is being experimentally tested in 
experiments where heavy nuclei are collided at high energies to produce such 
deconfined state of matter, the quark-gluon plasma. It is therefore imper- 
ative from the theoretical point of view to develop a better understanding 
of the behavior of gauge theories at temperatures above the deconfinement 
transition (for a recent review see e.g. Ref. PQ). 

Since a non-zero temperature for a quantum field theory is equivalent to 
a finite size and periodic boundary conditions in Euclidean time, in many 
important respects gauge theories at high temperature can be treated as ef- 
fectively living in one dimension less, a phenomenon generally referred to as 
dimensional reduction. More precisely, dimensional reduction can be defined 
through a specific pattern of pairwise degeneracies in the spectrum, char- 
acteristic of theories living in one less dimension than the original one, so 
that the occurrence of dimensional reduction can be quantitatively verified 
e.g. by Monte Carlo simulations. For £77(2) and SU(3) gauge theories in 
3 + 1 dimensions, this comparison was successfully performed in [2], where 
convincing evidence was provided that dimensional reduction occurs for tem- 
peratures all the way down to the deconfinement one. 

Much less is known about the actual form of the dimensionally reduced 
theory, especially for temperatures not too far from T c , where perturbative 
methods are not reliable. The simplest hypothesis, that the effective theory 
is the pure gauge theory with the same gauge group in one less dimension 
and in the confined phase, is ruled out by numerical data on the spectrum of 
screening masses for SU(2) [3 , while for SU(3) the agreement is better. A 
more sophisticated proposal was made in [3], where the effective 3D theory 
is taken to be a three-dimensional gauge theory coupled to an adjoint scalar. 
This effective theory has been shown to correctly reproduce the low-lying 
states of the spectrum for both ££7(2) [3] and SU(3) [S| for temperatures 
down to ~ 2T C . 

On the other hand, in the temperature range just above T c , we have the 
possibility of making a solid conjecture about the form of the dimensionally 
reduced theory, at least for the case of SU(2) or, in general, for gauge groups 
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such that the transition is of second order: this is the scaling region above 
the transition, where universality arguments hold and the effective theory 
can be conjectured to belong to the universality class of theories whose global 
symmetry group coincides with the center of the gauge group. For SU(2), the 
relevant universality class is the one of the 3D Ising model. This application 
of universality arguments to the deconfinement transition was introduced 
by Svetitsky and Yaffe 7^ and many of its predictions have been accurately 
verified: the values of the critical indices (the most recent results can be 
found in Refs. [3011011111), universal amplitude ratios [12], and correlation 
functions at criticality IH] . In this approach, the fundamental degree 
of freedom of the effective theory is the Polyakov loop, which is the order 
parameter of the deconfinement transition. 

It is far from obvious a priori that the higher masses in the spectrum 
should be universal, since only the lowest one contributes to the free energy. 
However, it has recently been shown that for the three dimensional Ising 
universality class, in the broken symmetry phase, there is strong evidence 
for such universality: the Ising model and lattice regularized 4 theory both 
exhibit a rich spectrum of massive excitations, and the mass ratios coincide 
in the scaling region [13] . Particularly interesting is the existence of an ex- 
cited mass in the scalar sector at less than twice the mass of the fundamental 
excitation: being the interaction attractive in the broken symmetry phase, it 
is natural to interpret such state as a bound state of two fundamental par- 
ticles, an interpretation that was confirmed by studying the Bethe-Salpeter 
equation for 3D </> 4 theory in the continuum |16j . 

It is therefore natural to ask whether the same spectrum characterizes 
SU(2) gauge theory in the scaling region above the deconfinement temper- 
ature. This paper is devoted to a numerical verification of this ansatz. In 
particular we will focus on the scalar sector: we will show that a bound 
state of the fundamental screening mass can indeed be extracted from the 
correlation functions of suitably defined operators, that are constructed from 
Polyakov loops in a manner similar to that employed in jT3] to study the 
spectrum of the Ising model. We also evaluated the lowest screening mass 
with angular momentum 2 in the scaling region and found a value compat- 
ible with the one of the Ising model. Preliminary results were reported in 
Ref. [13. 

In the following section we describe in detail the method we used and the 
numerical results, while Section 3 is devoted to the discussion of the results 
and some speculation about the case of SU(3). 
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2 Numerical determination of the spectrum 
of scalar screening masses 

As a first step we considered the connected wall-wall correlation function of 
the Polyakov loop on a lattice with N x x N y x N z x N t sites and lattice 
spacing a 

G(\ Zl -Z2\) = (P(zi)P(ft)> - <F(*i)><PteO> , 

where 



is the average of the Polyakov loop P(x,y } z) = Tr [] * =1 f/ 4 (x, y, z, 714a) 
over the xy-plane at a given z. The wall average implies the projection at 
zero momentum in the xy-plane. The correlation function G(\z\ — 22 1) takes 
contribution from all the screening masses in the + channel, i.e. from the 
lowest (fundamental) mass mi, from the first excited mass m.2, and so on. 
For a periodic lattice, we have 



G{\ Zl -z 2 \) = Q,e- miL 
+ C a 



g-mi|zi— 22! _j_ g-mi(L z -|2;i-22|) 



(1) 



+ C 2 



e -m2]zi— «2| _|_ e ~mi{L z -\zi-Z'2\) 



+ ... 



with L 2 = aA^. Here the dots represent the contribution from higher mass 
excitations and from multi-particle cuts. It should be noticed that, since 
(P(z)) 7^ in the high temperature phase, a ^-independent exponential term 
appears in the r.h.s. of Eq. (JH) [TH]. It is easy to verify that such term is 
sub-dominant with respect to the contribution of the fundamental mass and 
of the excited masses mj, % > 1 for which m^ < 2mi. 

We discretized the theory on the lattice with the Wilson action and 
generated Monte Carlo configurations by the overheat-bath updating algo- 
rithm [20] with Kennedy- Pendleton improvement [21]. Measurements were 
taken every 10 upgradings and error analysis was performed by the jackknife 
method applied at different blocking levels. 

We made simulations at two values of (3 slightly above the critical value 
P c = 2.29895(10) determined in Ref. 02|: at (3 = 2.33 (corresponding to 
T/T c ~ 1.10) we collected 1M measurements on a 18 2 x 36 x 4 lattice and at 
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(3 = 2.36 (corresponding to T/T c ~ 1.21) we collected 1.5M measurements 
on a 18 3 x 4 lattice. 

To make sure that our Monte Carlo ensembles were not contaminated by 
configurations with mixed phase, we plotted histograms of the Polyakov loop 
averaged over the lattice configuration. In all cases, histograms showed two 
non-overlapping peaks located on opposite sides with respect to the origin, 
thus indicating that no configurations with mixed phase were present in the 
thermal equilibrium ensemble. 

We determined the effective mass through 



G(z) = C e - m ^ Lz + C e~ m ^ z + e - m ^- z) (2) 

and found for both (3 values that data for function of z reach a 

plateau value, corresponding to the fundamental mass mi, as z is increased 
towards N z /2. The deviations from the plateau value at small z can be 
attributed to lattice artifacts and to the possible effect of other physical 
states. In order to single out the latter contribution, we rescaled the values 
of Tn ef f and z by mi = l/£ and respectively, and put together on the same 
plot data from f3 = 2.33 and (3 = 2.36. As shown in Fig.Q data from the two 
different /3's fall on the same curve before reaching the plateau, indicating 
that the pre-asymptotic behavior of the effective mass is not a lattice artifact, 
but it is instead a physical effect, due to higher mass excitations. 

Universality arguments hold in the scaling region of the deconfinement 
transition: in order to check whether our choice of (3 is in fact within this 
region, one can compare the values of the lowest screening mass mi at f3\ = 
2.33 and (3 2 = 2.36, whose ratio is predicted by universality to scale as: 

wni(ft) ( (3\- (3 C x ^ 



ami(/3 2 ) \(3 2 -(3, 

where v ~ 0.63 is the correlation length critical index for the three dimen- 
sional Ising model. The values we measured 1 give 

0^^ 0.3654(31) 

amxiPi) 0.5431(61) y ' y ' 



lr rhe value quoted for (3\, which is the point where the analysis of the whole spectrum 
was performed, was obtained with the variational method described below. The value at 
02 was obtained simply from the long-distance behavior of m ef f as defined in Eq. J2J. The 
value obtained with the latter method at (3\ is compatible with the one quoted. 
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Figure 1: Effective screening masses in physical units vs the separation be- 
tween walls in physical units at (3 — 2.33 on a 18 2 x 36 x 4 lattice and (3 = 2.36 
on a 18 3 x 4 lattice. 



to be compared to 

'A - A 



We conclude that the scaling region includes, within the accuracy of our 
data, these two values of /3, and hence that such region extends at least 
up to temperatures ~ 1.2 T c . On the other hand, the value of the mass at 
T = 2T C reported in Ref. [3] is rather far from the scaling prediction, therefore 
the scaling region does not extend up to 2T C . 

In order to detect possible excited states, we adopted the so-called "vari- 
ational" method [221 12H]- The method consists of the following steps: (i) 
defining a basis of (wall-averaged) operators {O a }, (ii) computing the con- 
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nected cross-correlation matrix among these operators, 



C aP {z) = (O a {z)O p {0)) - (O a (z))(Of,{0)) , (6) 

(iii) diagonalizing C a p{z) to obtain the eigenvalues \{z) and finally (iv) 
extracting the masses rrii through the relation 

X^z) = c + c \e~ miZ + e - m ^- z) \ . (7) 

The fundamental mass mi corresponds to the leading eigenvalue Ai, the 
first excited mass m 2 corresponds to the next-to-leading eigenvalue A 2 , and 
so on. In practical simulations, for each value of z we have a numerical 
determination of an "effective" mass rrii(z) through Eq. (|7ji. As z increases, 
rrii(z) reaches a plateau value which we identify with the screening mass m^. 
In our simulations, the signal-to-noise ratio allowed to extract a screening 
mass only from the two leading eigenvalues. 

The effectiveness of this method relies on the choice of a "good" set of 
operators; it is convenient to define operators living on different length scales 
(for instance, by use of recursive definitions). Moreover, by defining operators 
with non-trivial transformation under spatial rotation, it is possible to look 
for states with non-zero angular momentum. 

As a first set of operators with + quantum numbers, we considered those 
which can be built by adapting to the present case the recursive "smooth- 
ing" procedure used in Ref. [T3] for the 3D Ising model and the 3D lattice 
regularized 4 theory: 

P {0) (x,y,z) = P(x,y,z) , P {n+1) (x, y, z) = sign(w) (1 - w)\u\ + wy , (8) 

where u is the average of four nearest values of the Polyakov loop on the 
:q/-plane, 

u = - (p( n \ x -a, y, z) + P (n \x, y-a, z) + P (n \x+a, y, z) + P {n \x, y+a, z) 

and we chose w = 0.1 and y = (P(x,y, z)) . We considered a subset of six 
such (wall-averaged) operators corresponding to n = 0, 3, 6, 9, 12 and 15 
smoothing steps. At = 2.33 (T/T c ~ 1.10) on a 18 2 x 36 x 4 lattice, we 
collected ~1.6M measurements. 
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The relevant observable to check our ansatz is the ratio between the first 
excited and the fundamental screening masses in the scalar channel, which 
in the 3D Ising and <ft A theories assumes the value 

— = 1.83(3) . (9) 
mi 

The effective values of this ratio for SU(2) gauge theory are shown in Fig. 2: 
while the determination of this ratio is much more difficult than in the Ising 
model, the compatibility of the result with the expectations from universality 
is satisfactory. We have also observed that the excited scalar state cannot be 
detected when only the first four operators of our set (0, 3, 6 and 9 smoothing 
levels) are included in the analysis. This fact suggests that the first excited 
state in the + channel couples with operators living at rather large length 
scales. 
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Figure 2: Ratio of the first excited effective screening mass in the + channel 
to the fundamental one vs z at (3 = 2.33 on a 18 2 x 36 x 4 lattice. The 
statistics of the simulation is ~1.6M. 

Finally, in order to look for states with non-zero angular momentum, we 
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considered a second set of 2 + operators, inspired by Ref. [T8] . 

N X Ny 



P""(:) , P(nio, n 2 a,z) 



x 



w ni=l ri2=l 

P(n\a + na, ra 2 a, 2;) — P(nia, n 2 a + na, 2) , (10) 



with n = 1, . . . , 5. By the same procedure adopted in the case of the screening 
mass in the + channel, we got for the fundamental mass in the 2 + channel the 
value m!f ■* = 0.86(21)a _1 . Again, we studied the observable m^ 2 ^ /mf \ 
i.e. the ratio between the fundamental mass in the 2 + channel, determined 
from this second basis of operators, and the fundamental mass in the + 
channel, determined from the first basis of operators. The result of our 
determination, as a function of z, is shown in Fig. ©. Despite the noisiness, 
one can again see that the result is compatible with the value found [TH] in 
the 3D Ising universality class 2 

(2+) 

1 =2.56(4). (11) 

777.1 



3 Discussion 

The numerical results we have obtained are in agreement with our starting 
hypothesis: the spectrum of screening masses for SU{2) gauge theory near 
T c is indeed the one predicted by universality arguments, i.e. it coincides 
with the spectrum of the transfer matrix of the 3D Ising model in the scaling 
region of the broken-symmetry phase. In particular a bound state of two 
elementary screening masses can be detected. An unambiguous verification of 
the coincidence of the two spectra will require both the analysis of more states 
and a substantial increase in the accuracy in the masses we have determined. 

At first sight, application of these methods to the case of SU(3) does 
not seem possible, given that the deconfinement transition is first order. 
However, the transition is rather weak, both in the gauge theory and in the 3- 
state Potts model that would be in the same universality class if universality 

2 This mass ratio was actually evaluated in Ref. 18 for the 3D Ising gauge model: 
however duality implies that in the broken symmetry phase the spectrum coincides 
with the one of the 3D Ising model. 
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Figure 3: Ratio of the fundamental effective screening mass in the 2 + channel 
to the fundamental one in the + channel vs z at (3 = 2.33 on a 18 2 x 
36 x 4 lattice. The statistics of the simulation is ~1.8M. The horizontal 
lines represent the same mass ratio for the 3D Ising model |18| . with its 
uncertainty. 

arguments applied. Therefore the three-dimensional correlation length at 
criticality, while not diverging, assumes rather high values: hence, there 
exists a range of temperatures in which the correlation length is large in 
terms of lattice spacing, and it is not unreasonable to expect the spectra of 
the two theories to be similar in such region. We plan to investigate this issue 
studying the spectrum of the three-dimensional 3-state Potts model in the 
broken symmetry phase, and comparing it to the spectrum of SU(3) pure 
gauge theory near T c . 

Acknowledgement. We are grateful to M. Caselle for many useful discus- 
sions. 



10 



References 

J. B. Kogut, D. Toublan and D. K. Sinclair, |arXiv:hep-lat/0208076| 



S. Datta and S. Gupta, Nucl. Phys. B 534 (1998) 392 [arXiv:hep-lat/98- 
06034]. 

S. Datta and S. Gupta, Phys. Lett. B 471 (2000) 382 [arXiv:hep-lat/99- 
06023]. 

K. Kajantie, M. Laine, K. Rummukainen and M. E. Shaposhnikov, Nucl. 
Phys. B 503 (1997) 357 |arXiv:hep-ph/9704416| . 

A. Hart and Q. P hilipsen, Nucl. Phys. B 572 (2000) 243 
|arXiv:hep-la t /9908041] . 

A. Hart, M. Laine and O. Philipsen, Nucl. Phys. B 586 (2000) 443 



[2] 
[3] 
[4] 
[5] 
[6] 

|arXiv:hep-ph70 004060 . 

[7] B. Svetitsky and L. G. Yaffe, Nucl. Phys. B 210 (1982) 423. 

[8] S. Fortunato, F. Karsch, P. Petreczky and H. Satz, Nucl. Phys. Proc. 
Suppl. 94 (2001) 398 |arXiv:hep-lat/0010026| . 

[9] R. Fiore, A. Papa and P. Provero, Nucl. Phys. Proc. Suppl. 106 (2002) 
486 |arXiv:hep-lat/0110017 . 



[10] R. Fiore, A. Papa an d P. Provero, Phys. Rev. D 63 (2001) 117503 
|arXiv:hep-lat/0102004|. 



[11] A. Papa and C. Vena, |arXiv:hep-lat/0203007| 

[12] J. Engels and T. Scheideler, Nucl. Phys. B 539 (1999) 557 
|arXiv:hep-lat /9808057| . 

[13] F. Gliozzi and P. Provero, Phys. Rev. D 56 (1997) 1131 
|arXiv:hep-lat/97010T4| . 

[14] R. Fiore, F. Gliozzi a nd P. Provero, Phys. Rev. D 58 (1998) 114502 
|arXiv:hep-lat/9806017| . 

[15] M. Caselle, M. Hasenbus ch and P. Provero, Nucl. Phys. B 556 (1999) 
575 |arXiv:hep-lat/9903011| . 

11 



[16] M. Caselle, M. Hasenbusch, P. Provero and K. Zarembo, Nucl. Phys. B 
623 (2002) 474 |arXiv:hep-th/0103130|. 



[17; 

[18 

[19 
[20 
[21 
[22 
[23 



R. Fiore, A. Papa and P. Provero, j|arXiv:hep-lat/0208020| 

V. Agostini, G. Carlino, M. Caselle and M. Hasenbusch, Nucl. Phys. B 
484 (1997) 331 |arXiv:hep-lat/9607029| . 

I. Montvay and P. Weisz, Nucl. Phys. B 290 (1987) 327. 

R. Petronzio and E. Vicari, Phys. Lett. B 254 (1991) 444. 

A. D. Kennedy and B. J. Pendleton, Phys. Lett. B 156 (1985) 393. 

A. S. Kronfeld, Nucl. Phys. Proc. Suppl. 17 (1990) 313. 

M. Luscher and U. Wolff, Nucl. Phys. B 339 (1990) 222. 



12 



